Topological Phases for Fermionic Cold Atoms on the Lieb Lattice 
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We investigate the properties of the Lieb lattice, i.e a face-centered square lattice, subjected 
to external gauge fields. We show that an AbeHan gauge field leads to a pecuhar quantum Hall 
effect, which is a consequence of the single Dirac cone and the fiat band characterizing the energy 
spectrum. Then we explore the effects of an intrinsic spin-orbit term - a non-Abelian gauge field - 
and demonstrate the occurrence of the quantum spin Hall effect in this model. Besides, we obtain 
the relativistic Hamiltonian describing the Lieb lattice at low energy and derive the Landau levels 
in the presence of external Abelian and non-Abelian gauge fields. Finally, we describe concrete 
schemes for realizing these gauge fields with cold fermionic atoms trapped in an optical Lieb lattice. 
In particular, we provide a very efficient method to reproduce the intrinsic (Kane-Mele) spin-orbit 
term with assisted-tunneling schemes. Consequently, our model could be implemented in order to 
produce a variety of topological states with cold-atoms. 

PACS numbers: 37.10.Jk, 05.30. Fk, 71.70.Di 



I. INTRODUCTION 

During the last decades, topology has influenced many 
fields of physics through the renewed description of var- 
ious phenomena. In condensed-matter physics, topologi- 
cal invariants — known as Chern numbers — have played 
an important role in the description of the integer quan- 
tum Hall effect (IQHE) [1]. Here, the quantized Hall 
conductivity of a two-dimensional (2D) electron system is 
expressed as a sum cth = J2e„<Ef ^chiE„) of Chern 
numbers N^hiEn) that are integers associated with the 
energy band En [2, 3]. Here Rk is von Klitzing's constant 
and Ep denotes the Fermi energy assumed to lie inside 
a gap of the bulk energy spectrum. Furthermore, it has 
been proven that the sum of Chern numbers is expressing 
the number of gapless edge states located inside the bulk 
energy gaps. These edge states carry the current in the 
IQHE [4, 5]. 

The breaking of time-reversal symmetry (TRS) due to 
external magnetic fields plays a crucial role for the topo- 
logical interpretation of the IQHE [6]. Recently, the dis- 
covery of the so-called quantum spin-Hall effect (QSHE) 
has lightened a new path for the investigation of sys- 
tems where TRS is preserved [7-12]. The QSHE mani- 
fests itself in insulating systems that show a non-trivial 
Z2 index [7]. These so-called topological insulators are 
characterized by the presence of spin-filtered edge states 
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in the gaps of the bulk energy spectrum. Because of 
TRS invariance, the spin-up and spin-down states move 
in opposite directions along the edge of the system. As 
a consequence, the total charge current as well as the 
associated Chern numbers are zero [12, 13]. Yet, a spin- 
Chern number has been introduced in order to measure 
the spin-transport [14] and to distinguish the Z2 class of 
the system [15]. 

The interplay between the lattice topology and the 
QSHE has been the focus of various recent investiga- 
tions. In particular, the QSHE has been studied for 
the Kagome [16, 17], the Lieb and Perovskite [18], the 
honeycomb [7, 19], the square [20, 21], the T3 [22], the 
checkerboard [23], the pyrochlore [24], diamond [25] and 
the square-octogon [26] lattices. 

In this context, some lattice models are of particular 
interest as they show dispersionless energy bands. These 
flat bands correspond to a macroscopic number of de- 
generate localized states. Originally, flat bands played a 
fundamental role in magnetism, as they were shown to 
accompany the occurance of ferromagnetic ground states 
in multi-band Hubbard models [27-29]. More recently, 
the existence and the robustness of these special bands 
have been extensively studied in a vast family of frus- 
trated hopping models [30, 31] and for the case of elec- 
tron localization due to magnetic fields and spin-orbit in- 
teractions [32] . Interestingly, singular touchings between 
fiat and dispersive bands have been shown to be topo- 
logically protected by real-space loops [30]. On the face- 
centered square lattice — also known as Lieb lattice — 
a fiat band touches two linearly dispersing bands, i.e.. 
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the flat band intersects a single Dirac point, and the low- 
energy regime is described by a quasi-relativistic equation 
for spin-1 fermions [33, 34]. 

Nowadays, various lattices can be engineered using 
cold atoms trapped by electromagnetic fields [34-39]. 
In particular, the realization of topological states of 
matter with cold fermionic atoms appears to be a 
realistic and attractive goal from the experimental point 
of view [20, 21, 40]. A significant advantage of these 
experiments is the full control of a wide range of system 
parameters as, e.g., lattice geometry, interaction and 
disorder. In these experiments, engineered gauge fields 
allow to mimic the effects of magnetic fields [41-43] or 
spin-orbit interactions (SOIs) [20, 21, 44-48]. These 
gauge fields can be generated by spatially-varying laser 
or magnetic fields which modify particle-hopping via 
non-trivial Berry's phases [49, 50]. Recent experiments 
have implemented light-induced external gauge fields 
and reproduced the physics of charges subjected to 
electric or magnetic fields [41, 42, 51]. Moreover, with 
such a setup one expects to observe several fundamental 
phenomena including the Hofstadter butterfly [43, 52], 
atomic analogues of the quantum Hall effects [53, 54], 
relativistic physics [55, 56], and vortex structures 
[41, 57, 58]. Optical-lattice setups also allow to consider 
a generalization of the ongoing experiments, namely the 
implementation of non-Abelian gauge fields [46, 59-61]. 
In particular, non-Abelian gauge fields acting on multi- 
level atomic systems could mimic SOI [44-46, 48, 62], 
paving the way to study the spin Hall [63] and quantum 
spin Hall effects [17, 21]. Very recently, a concrete 
proposal of an optical Lieb lattice for cold atoms has 
been presented [37]. In the later work, Apaja et al. have 
shown that a fermionic cloud expanding after the release 
of the harmonic trap should show clear signatures of 
the ffat band's localized states. Finally, the existence 
of flat bands with non-trivial topological order has 
been demonstrated [64], contradicting the belief that 
non-dispersive bands were associated to vanishing Chern 
numbers [31]. 

Motivated by the possibility to engineer an optical Lieb 
lattice for cold fermionic atoms, we investigate the emer- 
gence of topological properties for various conflgurations 
of synthesized gauge flelds. We first provide an original 
analysis of a peculiar IQHE, in the case where a uni- 
form magnetic field is present in the Lieb lattice. We 
then explore the effects of an intrinsic spin-orbit term [7] 
and show how it leads to quantum spin Hall states. In 
this framework, we extend the seminal work of Ref. [18] 
and derive an effective Hamiltonian describing the low- 
energy regime. This Weyl-like Hamiltonian leads to a 
three-component quantum equation that resembles the 
relativistic equation for spin-1 particles. Besides, we ob- 
tain the Landau levels in the presence of an external mag- 
netic field and spin-orbit interaction. Finally, we discuss 
the optical-lattice realization of this Lieb system and pro- 
pose realistic methods for creating Abelian (magnetic) 
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FIG. 1: The face-centered square lattice or Lieb lattice. The 
Peierls phases e'*'™-', where 6{m) = n^m, are associated to a 
uniform magnetic flux per plaquette $ and are indicated by 
vectical black arrows. We set x = 2m£o and y = 2n£o. 

and non-Abelian (spin-orbit) gauge fields. We show that 
the Lieb lattice is particularly suited to reproduce the in- 
trinsic spin-orbit term introduced by Kane and Mele [7]. 
The later, which involves complex spin-dependent nest- 
nearest-neighbour hoppings, can be simply decomposed 
into nearest-neighbour hopping on a square sublattice. 
This elegant idea is a non-Abelian generalization of the 
method proposed in Ref. [17] for generating the Abelian 
Haldane-type gauge field. 

II. THE LIEB LATTICE AND TOPOLOGICAL 
PHASES IN EXTERNAL FIELDS 

We consider the face-centered square (Lieb) lattice, 
which is shown in Fig. 1. This lattice has a unit cell char- 
acterized by three lattice sites, hereafter referred to as H, 
A and B. Site H has four nearest-neighbors (NN), namely 
two A and two B sites. On the contrary the A and B sites 
have only two NN H sites. The bulk properties of the Lieb 
lattice can be analyzed within a tight-binding (TB) ap- 
proximation. In this limit, the Hamiltonian of the system 
can be written as Tio = ^ j)a '^L'-ja with spin inde- 
pendent NN hopping amplitude t. Here, cj^(cjQ) is the 
creation (annihilation) operator for a particle with spin 
direction a on the lattice site j. In absence of external 
fields the problem can be diagonalized exactly and the 
spectrum reads 

eo(fc) = 0, (la) 
e±{k) = ±iV4 + 2cos(i;i • k) + 2cos(i;2 • k), (lb) 

where k — {kj;,ky) and are the lattice vectors, c.f. 
Fig. 1. The bulk energy spectrum is shown in Fig. 2a — it 
depicts two identical, electron-hole symmetric branches 
£-1-. Moreover, the Lieb lattice presents a unique non- 
dispersive band at the charge neutrality point (CNP). 
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FIG. 2: Energy spectrum of the Lieb lattice as a function 
of the momentum k = (k^, ky). The dashed lines delimit the 
first Brillouin zone. Panel (a): spectrum without external 
fields. Panel (b): energy spectrum in the presence of spin- 
orbit interaction tso = 0.1 t. 

This band is rooted in the lattice topology, which allows 
for insulating states with finite wave function amplitudes 
on the A and B sites and vanishing amplitudes on H sites. 
This property holds also when hopping to higher-order 
neighbors is allowed. Note that the three bands touch 
at the center of the first Brillouin zone, which we set for 
simplicity at F = 7r/2£o(l,l). The resulting properties 
of carriers in proximity of the F point are investigated in 
Sec. in. 



A. Uniform magnetic field and quantum Hall 
phases 

We now study the effects of a uniform magnetic field 
B = Bz on the spectral and transport properties of the 
Lieb lattice. We consider the Landau gauge 

A=(0,Sx,0)= fo,— ,0j , (2) 

where $ = 4>^^ B ■ AS is the number of magnetic flux 
quanta $o per plaquette and x = 2m£o. Hereafter we 
use the notation {m,n,(), with ( ~ {A,B,H}, to label 
the lattice sites. The gauge field A modifies the hopping 
along the y direction through x-dependent Peierls phases 
t — ^ ie'^*^™-', where the phase reads 

p{m,n-\-l.Ii) /'(m.n.B) 

e{m)^ / A-dl^ / A-dZ = 7r$m, (3) 

J(m,n,H) J{rn,nM) 

as illustrated in Fig. 1. Here, the integrations are per- 
formed along the links connecting the neighboring B and 
H sites. 

Setting $ = p/q, where p and q are mutually prime 
integers, the system becomes g-periodic along the x di- 
rection. By considering periodic boundary conditions. 



it is possible to diagonalize the resulting 3g x 3q spec- 
tral problem. This leads to the fractal energy spectrum 
shown in Fig. 3. As a function of the flux 3", the allowed 
energies depict two Hofstadter butterflies separated by a 
flat band at = [65] . This specific band is reminiscent 
of the flat band obtained at zero magnetic field. 

The fractal energy spectrum of lattices subjected to 
uniform magnetic fields are intimately related to the 
IQHE [67, 68]. When the Fermi energy Ep is located 
in a spectral gap, the Hall transverse conductivity of 
the system is quantized. This relation is supported by 
a Diophantine equation [69] which expresses the quan- 
tized Hall conductivity in terms of the magnetic flux and 
the position of the gap: In the r-th gap, the Hall con- 
ductivity is given by a^y — (e'^ /h)tj., where the integers 
{tr, Sr) satisfy 

r—ptr + qsr- (4) 

In general, the solutions (t^, Sr) are not unique and addi- 
tional criteria are needed in order to find quantized values 
of the Hall conductivity [69]. The integer also has a 
topological interpretation, since it represents the sum of 
Chern numbers characterizing the bands below the Fermi 
energy Ep: 

U^- J2 Nch{Ex). (5) 

Ex<Ef 

In order to investigate the quantum Hall phases in the 
Lieb lattice, we have numerically evaluated the Chern 
numbers Nch{E\) using the method of Ref. [70]. We 
have verified that the integer tr satisfies the Diophantine 
equation with the specific condition \tr\ < q/2 [69]. The 
full phase diagram describing the integer quantum Hall 
effect for the Lieb lattice is drawn in Fig. 3. It represents 
the infinitely many quantum Hall phases, characterized 
by the quantized transport coefficient cr^y — (e^//i)tr, 
inside the spectral gaps. The different positive [resp. 
negative] values of the Hall conductivity are designated 
by cold [resp. warm] colors. 

To identify the Hall plateaus stemming from the uni- 
form magnetic field, we represent the Hall conductivity 
f^xyiEp) as a function of the Fermi energy in the low-fiux 
regime $ ^ 1 (cf. Fig. 4). In this regime, the quantized 
conductivity evolves monotonically but suddenly changes 
sign around the van Hove singularities (VHS) [71], lo- 
cated a.t E = ±2 t (see the alternation of cold and warm 
colors in Figs. 3 and 4). Note that the gaps surrounding 
the topological fiat band at i? = 0, correspond to normal 
band insulators with vanishing conductivity axy = 0, for 
all values of the fiux $. This is a consequence of the fiat 
band's vanishing Chern number [31]. Most importantly, 
we observe that the Hall sequence presented in Fig. 4 
shows steps of Aa^y — (e^/h). It is interesting to com- 
pare the latter result with the Hall sequences obtained 
for the Ta and honeycomb lattices [22, 71], which are 
characterized by steps of Aaxy = 2(e^//i) between the 
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FIG. 3: Spectrum E — -E($) and phase diagram for $ — p/q 
with q < 47. The eigenvalues are dark blue dots forming two 
successive butterflies. Gaps are filled with cold [resp. warm] 
colors according to the related positive [resp. negative] values 
of the quantized conductivity a^y. The white gaps located 
around E = correspond to a^y = 0. 
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FIG. 4: Hall conductivity a^yiEp) as a function of the Fermi 
energy for "I> = 1/11. Cold [resp. warm] colors correspond to 
positive [resp. negative] values of the quantized conductivity, 
while the central plateau at a^y = is represented in magenta. 



VHS. This major difference [72] is due to the fact that 
the Lieb lattice is characterized by a single Dirac-Weyl 
point, whereas the Ta and honeycomb lattices display two 
Dirac cones. 



B. Spin-orbit interaction and quantum spin Hall 
phases 

In this Section, we study the effects of spin-orbit in- 
teractions (SOIs) on the properties of the Lieb lattice. 
Specifically, we introduce an intrinsic SOI term in the 
TB Hamiltonian, in analogy with the model of Kane and 
Mele for graphene [6-8, 18, 22, 73]. This term is mod- 
eled via a spin-dependent next-nearest-neighbor (NNN) 
hopping term 

"Hso = i^so XI X! "^U ^ ■ ^^-l^ ■ 

The (Tap are matrix elements of the Pauli-matrices cr 
with respect to the final and initial spin states a and 
/? and di/j are the two displacement vectors of the NNN 
hopping process connecting sites k and I. Since in 2D lat- 
tices hopping is naturally restricted to in-plane processes, 
the SOI is effectively proportional to cr^. Because of the 
unequal connectivity of A/B and H sites, the term (6) 
effectively induces hopping between A and B sites only, 
i.e., other NNN-hopping processes cancel. The spectrum 
is obtained by exact diagonalization [18] and reads 

er\k)^0 (7a) 

e^i°\k) = 2 [t^ (cos2(i;i • k) + cos2(t;2 • k)) (7b) 

+4:tlo sin^{vi ■ k) si-a^{v2 ■ k)]^ . (7c) 

The bulk energy spectrum is shown in Fig. 2b for tgo = 
0.1 1. Due to its topological origin, the non-dispersive 
band at £^ = is not affected by SOI. However, this 
term opens two bulk energy gaps Agap — 4 ^so between 
the non-dispersive and the electron/hole branches, re- 
spectively. 

The SOI has dramatic consequences on the transport 
properties of the Lieb lattice: as shown by Weeks and 
Franz, the gap Agap allows for a topological insulating 
phase [18]. The latter is characterized by a robust spin 
transport along the edges of the system. This quantum 
spin Hall phase, induced by the Haldane-type term (C), 
can thus be directly visualized when studying a finite 
piece of the lattice, i.e., by considering its edges. 

A standard method consists in diagonalizing the TB 
Hamiltonian with periodic boundary conditions imposed 
along one of the spatial directions. This abstract cylinder 
contains two edges and already allows to demonstrate 
the existence of helical edge states induced by the SOI. 
The corresponding energy spectrum (c.f. Fig. 5a) depicts 
several edge-state channels: for each energy value within 
the bulk energy-gap, there exists a single time-reversed 
(or Kramers) pair of eigenstates localized on each edge of 
the lattice. The conservation of TRS prevents the mixing 
of this couple of states by small external perturbations 
and scattering from disorder [7, 8]. 

The helical edge states characterizing the QSH phase 
are topologically protected against external perturba- 
tions. Their property can be quantified by looking at 
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the Z2-index v [7]. This topological invariant character- 
izes the eigenstates defined in the bulk — it is defined on 
a 2-torus, in direct analogy with the Chern numbers in- 
troduced in the quantum Hall effect. Following Ref. [f8], 
we have calculated this Z2-index v using the inversion 
symmetry of the lattice [25]. We obtain that SOI opens 
a spectral gap characterized by the index v — \^ there- 
fore classifying the Lieb lattice as a quantum spin-Hall 
insulator. 

In the absence of spin-mixing perturbations, the Z2 
index is related to the spin Chern number n„ [14, 15], 
through the simple relation v — nCTmod2, where ric — 
{N^ — Ni)/2 and N^^i represent the Chern numbers as- 
sociated to the individual spins. Using the numerical 
method of Ref. [70], we have obtained n^j — 1 in agree- 
ment with the above result. 

It is interesting to extend the analysis above by con- 
sidering the more realistic open geometry, i.e. a finite 
piece of Lieb lattice, thus characterized by a unique 
edge. We have solved the TB problem for a Lieb lat- 
tice of 39 X 39 sites with realistic straight edges. This 
yields a discrete energy spectrum extending in the range 
E e [-2.8t,2.8t]. In the vicinity of the CNP, within 
a range corresponding to Agap = 4 ^so, the eigenvalues 
correspond to eigenstates that are localized at the edge of 
the system. This result is illustrated in Fig. 5b, where the 
amplitude \'ijj-f{x,y)\'^ is drawn for a particular edge-state 
a.t E = 0.5t. Note that this coincides with \tpi{x,y)\'^ . 
Using this geometry, one can verify the helical property 
of these edge-states by computing their associated cur- 
rent, which for spins cr can be expressed as 

j^(m,n) = -iio^ip*{m + l,n)Ur,il;„{m,n) (8) 
+ V2ilj*{m + l,n + l)Vi;„{m,n) - h.c^ , 



(9) 



+ V2ip* {■m + l,n + l)Vip^{m,n) - h.cj . 



Here ?/'o-(TO,n) = {tlja-(m,n,ll),'ip^{rn,n, A),^pa■{rn,n,B)) 
and 
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We have verified that the currents J^(to, n) and j^{m, n) 
associated to the edge-states yield two vector fields cir- 
culating along the edge in opposite directions. 




FIG. 5: Panel (a): Energy spectrum E — E{ky) in the cylin- 
der geometry with tso = 0.2t. The bulk gap is traversed by 
gapless edge-states [i.e. helical states] . Panel (b): Edge-states 
amplitude \jp^{x,y)\^ for the open Lieb lattice with 39 x 39 
sites, straight edges and tso = 0-5t. This localized eigenstate 
satisfies Hxf}x{x,y) = Extp\{x,y) and corresponds to the en- 
ergy E\ — 0.5t which Hes within the gap. 



III. LOW-ENERGY FERMIONS: THE 
QUASIRELATIVISTIC REGIME 

In this Section, we focus on the the properties of the 
Lieb lattice for non-interacting fermions at low energy, 
i.e. close to the CNP. We perform a long wavelength ap- 
proximation to the Schrodinger equation underlying the 
TB Hamiltonian which consists in expressing the spatial 
part of the wave function as the product of a fast-varying 
part times a slow- varying part. Within this approxima- 
tion, the wave function can be written as 



(11) 



where a G {A,B,H} and Ra is the lattice site coordinate. 
We substitute this wave function into the Schrodinger 
equation and expand the slow-varying part as 

F^{R^,±dj) ~ F^{R^,)±dyVr F^{r)l^^ , +0(1^1') ■ 



Collecting all terms, we are left with 



T-L = -wpS • p , 



(12) 



where v-p = 2£ot 
iPx,Py,0) = -iti{dx 
trices are defined as 



is the Fermi velocity, and p — 
, 9y,0). Here the pseudo-spin ma- 






(13) 



These matrices fulfill the algebra of the angular momen- 
tum [EijSj] = icy-fcEfc and form a 3-dimensional repre- 
sentation of SU(2). However, contrary to the Pauli matri- 
ces, they do not form a Clifford algebra, i.e., {S^, Ej} ^ 



25. 



Therefore, while Eq. (12) describes electrons with 
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a linear energy spectrum, it does not represent a Dirac 
Hamiltonian. By introducing a rotation operator around 
the z axis defined by Vzicj)) = exp (— iE^i/)), a generic 
state |a) is transformed into itself by 'Dz{2n)\a) -> |a), 
implying that the pseudo-spin S describes an integer spin 
5=1. 



A. Spin-orbit interaction 

Within the long wavelength approximation we can also 
express the intrinsic SOI introduced in the previous Sec- 
tion. This term reads 



(14) 



where Ago is the effective spin-orbit coupling strength 
and is a Pauli matrix. The energy spectrum can be 
computed in this regime and reads 



;.(SO) 



= 



?(so) _ 



so ■ 



(15a) 
(15b) 



This is two-fold degenerate, with degeneracy correspond- 
ing to spin-up and spin-down. 



with corresponding eigenvalues e±^n = ±A/2n + 1 and n 
integer. Here the 0n for n > are the eigenfunctions of 
the one-dimensional harmonic oscillator, 



0n(e) 



1 



V2"7ri/2n! 



-e/2 



(19) 



where hn denotes the Hermite polynomial of order n, 
while we define 0_i = 0. In addition to the eigenfunc- 
tions (18) there are two different types of solutions at 
energy £o,n = 0. The first of these is related to the flat 
band at zero magnetic field (la) and reads 



1 



V'0,n(O = ^ I I , (20) 



where n > 0. The second type of zero-energy Landau 
level is given by 




(21) 



Note that tpl^\^) is not a generalization of (20) to the 
case n = 0. 



C. Spectrum with magnetic field and spin-orbit 
interaction 



B. Landau levels in a uniform magnetic field 

The long wavelength approximation also allows to com- 
pute the Landau levels that arise in the presence of the 
uniform magnetic field. The system Hamiltonian reads 



-Hb =wfS- [p- -A 



(16) 



where A is the vector potential associated with the mag- 
netic field B = (V X A) perpendicular to the lattice 
plane. 

We solve the Schrodinger equation in the Landau gauge 
with A = {—By ,0,0). Further, we make the Ansatz 
vj/ = '0(y) exp(ifcx). Introducing a /c-dependent shift in 
the ^-coordinate, \/B£, = By + k we are left to solve a 
system of coupled linear differential equations 



(17a) 
(17b) 
(17c) 



^Vh(0 = e^A(0 
eV'A(0-i^B(0=£~V'H(C) 



for the three components of '!/'(C)- Here e = E / [hv-pv B) 
is the rescaled eigenenergy. The Landau levels at non- 
zero energy are given by 



^±,n(e) 



i)/^/2 



(Vn0n-1 + \/n+ l(?!)n4 

iVStm^n I (18) 



Now we turn to the effects of finite SOI on the Landau 
levels obtained above. The energies e^.o-.n (with a = 
{0, ±} and (J = {t, i}) of the Landau levels are the three 
solutions of 



+ — 



Vn+T 



£«, cr, n — cr^SO £Q,(T,n + fAso 

The corresponding wave functions read 



(22) 




In the case of weak SOI, i.e. Ago ^ 1, the Landau 
levels are given by 

±V2^- ^^ + 0{Alo) (24) 



£±,(7,11 

£0,(T,n 



aA 



so 



2n+ 1 



4n + 2 
so) 



(25) 



Consequently, the main effect of finite SOI is to lift the 
spin degeneracy, with a level separation that decreases 
with growing Landau level index n. Moreover, the former 
highly degenerate zero energy levels, c.f. Eqs. (20) and 
(21), are now split into a family of fiat bands at energies 

£0.(T,n- 
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FIG. 6: (Color online) (a) Distribution of the laser field in- 
tensity for generating a Lieb lattice, (b) Two cuts through 
(a) corresponding to the B and H sites (red line) and A sites 
(blue line) 



IV. OPTICAL LATTICE REALIZATION 

Experimentally, the Lieb lattice can be realized as an 
optical lattice created by six counter propagating pairs 
of laser beams. Four pairs are aligned along the x and 
y directions, two with wavelength A = a-nd two with 
wavelength A = Finally, it requires two other laser 

pairs with a direction of ±45° with respect to the x- 
axis. A detailed procedure leading to this choice of laser 
configuration has been discussed in Refs. [34, 37]. The 
potential profile is given by the field 



VoL{x,y) =Vo [sm^{kx) + sin^{ky)] 

+ Vi [sin^(2fcx) + sin^(2fc?/)] 

r + y' 



(26) 



cos 



cos 



x-y 



with Vq — Vi = 2V2 and k = tt/A, which is depicted 
in Fig. 6. It is apparent that the hopping probability 
between A and B sites is exponentially small compared 
to the hopping probability between neighboring H and 
A/B pairs. 



A. Simulation of the U(l) synthetic gauge field and 
quantum Hall phases 

Recently, synthetic U(l) magnetic fields for cold neu- 
tral atoms have been proposed [43, 50, 74] and experi- 
mentally realized [75]. In such setups, atoms reproduce 
the dynamics of charged particles subjected to a uniform 
magnetic field and can effectively show quantum Hall 
phases. Several methods can be used in order to simu- 
late the Hofstadter model [52] with these systems. These 
methods are generally based on the fact that the Peierls 



phases can be engineered by external optical [43, 74] or 
magnetic [21] fields. These electromagnetic fields can in- 
deed induce hopping between neighboring lattice sites, 
when the latter host atoms in different internal states, 
say \g) and ]e). More precisely, the external fields trigger 
(Raman) couplings between these internal states, result- 
ing in a NN-hopping amplitude 



tr. 



(27) 



where the Rabi frequency ^g^e typically includes space- 
dependent phase factors and where we suppose the states 
\g) and ]e) to be trapped in neighboring sites [43, 74]. 
Here, the Peierls phase exp{i9{xg)) is directly related to 
the coupling laser's wave vector. 

The gauge field (2), which leads to the Peierls phase 
(3) can be readily engineered in an optical lattice ex- 
periment by exploiting these methods [43, 74]. From 
Fig. 1, it is clear that the a;-dependent assisted- hopping 
involves nearest-neighbors and occurs along the y direc- 
tion only, i.e. between B and H sites. In this sense, 
the phases 9{m) = ir^m can be realized on the opti- 
cal Lieb lattice by extending the methods envisaged for 
the standard square lattice (cf. Ref. [43, 50, 74]): In 
the Lieb lattice case, one should trap two internal states 
alternatively along the y-direction: the B-sites (resp. H- 
sites) should host an atom in the internal state \g) (resp. 
]e)). Coupling these states with external fields should 
then induce hoppings of the form (27), resulting in the 
space-dependent Peierls phase (3). Note that for gener- 
ating a magnetic flux $ per plaquette, a double phase 
d\j = 26 ^ 27r$m is required for the square lattice com- 
pared to the Lieb lattice. 

In the previous sections, we have discussed the exis- 
tence of quantum Hall phases in a fermionic Lieb lattice 
subjected to a uniform magnetic field. In order to pro- 
duce these phases in a cold atom experiment, one should 
engineer a U(l) gauge field for fermionic atoms. As al- 
ready discussed in Ref. [21], most of the schemes gener- 
ating gauge fields for bosons use Raman transitions to 
couple the internal states, and would lead to high spon- 
taneous emission rates for fermionic atoms. Therefore 
realizing (integer) quantum Hall states would require al- 
ternative methods [61, 76, 77]. Such a proposal was 
introduced in Ref. [21] and uses radio- frequency mag- 
netic fields produced by a set of current-carrying wires. 
The latter are periodically spaced on an atom chip and 
drive transitions between several internal states of ^Li 
fermionic atoms. These effective "Raman transitions" 
lead to assisted hopping (27) and can be tuned in order to 
produce the desired Peierls phases. In order to engineer 
the U(l) gauge field (2), one can simplify the method ini- 
tially proposed in Ref. [21] (which leads to the creation of 
SU(2) gauge fields) and consider transitions between two 
internal states of ^Li, e.g. \gi) = \F~l/2,mF = i/'2.) and 
]ei) — ]3/2,l/2). We stress that this practical scheme 
can be directly generalized to the Lieb lattice in order 
to generate the phases 9{'m) = 7r$m accompanying the 
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FIG. 7: The Lieb lattice with the intrinsic spin-orbit cou- 
phng. This coupHng is equivalent to a non-Abelian gauge 
field A leading to spin-dependent Peierls phases, as indicated 
by red arrows. 



hopping between neighboring H and B sites. 



B. Simulation of the SU(2) gauge field with 
neutral atoms 

Here, we describe a practical scheme to simulate the 
intrinsic SOI term in a fermionic Lieb lattice. This 
coupling is equivalent to a SU(2) gauge field and could 
therefore be engineered in a multi-component atomic 
system through state-dependent Peierls phases (cf. 
Fig. 7). Our proposal is based on the observation that 
the spin-dependent NNN-hoppings are equivalent to 
AW-hoppings defined on the square sublattice formed 
by the A and B sites only [40] (cf. Fig. 8 (a)). Con- 
sequently, generating the SOI reduces to the simple 
problem of engineering the Peierls phases e^^/'^Uz on 
a rotated square lattice, which we now label using the 
notations m and n (cf. Fig. 8 (b)). Obviously, the 
subtlety relies in the orientation of these phases: the 
phases are positive for a particle hopping respectively 
clockwise and anti-clockwise in neighboring plaquettes 
(cf. Fig. 8 (a)). Note that this important fact naturally 
leads to a staggered magnetic field, with fluxes ±$, for 
each spin component (cf. Fig. 8b). We note that, in 
order to reproduce such a staggered field, one can simply 
exploit the fact that the hopping induced by Raman 
transitions between the internal states g and e is such 
that (te*^(=^s))e,g = (te^®*=^''^);,e [61]. 

Let us now describe a feasible and concrete scheme to 
synthesize the SOI term (6) in a cold-atom experiment. 
Our proposal requires four states ei,2 and and ex- 
ternal fields producing Raman transitions in both the x 
and y directions (cf. Figs. 8 (a)-(d)). Such states can 
be chosen as being four internal states \F,rap) of ^Li, 
e.g. =11/2,1/2), I52) =13/2,-1/2), |ei) =[3/2,1/2), 



(a) (b) 




m-2 m-1 m m+1 m+2 n-2 n-1 n ni+2 



FIG. 8: (a) The Lieb lattice with the intrinsic spin-orbit 
coupling. The Peierls phases (red arrows) are defined on 
the square sublattice formed by the A and B sites, (b) 
Laser assisted- hoppings (arrows), associated couplings and 
synthetic fiuxes felt by the up-spin component. Four states 
ei,2 and g\^2 are trapped in state-dependent potentials and 
are represented by filled or empty orange or red disks (cf. fig- 
ures below), (c)-(d) State-dependent potentials and couplings 
along the x and y directions. 

and [62) =[1/2, —1/2) [21]. First, one needs to trap these 
states in state-dependent lattices [21, 59, 61] along the x 
and y directions (cf. Figs. 8c-d). Then external fields 
should drive Raman transitions between these states, 
with the corresponding Rabi frequencies 



X 

31 .ei 






— 0^ 


V 

ei,gi 




62, 32 


- 



as indicated by arrows in Figs. 8c-d. At this point, 
we emphasize that the Rabi frequencies are controlled 
by the coupling lasers and that they are chosen to 
be different for transitions driven along the x and y 
directions. Note that these Rabi frequencies typically 
contain phase factors depending on the coupling lasers 
wave vectors [43]. Moreover the Rabi frequencies 
associated to the opposite transitions are simply given 
by ^^^,g^ = i^^,,e,y^ where fi = i,y. Therefore, 
the hopping amplitudes along a given direction are 
accompanied by Peierls phases with alternating signs 
[61]. This leads to a staggered magnetic field for each 
spin component, with fluxes ±$, as illustrated in Fig. 8b. 

Now, the SOI term (6) requires a specific configuration 
of these Rabi frequencies. The associated SU(2) gauge 



field is proportional to ctz, which is simply achieved by 
imposing the constraint $7^ = (^^2)* (i-^- the coupling 
lasers should be characterized by opposite wave vectors) . 
Besides, the desired gauge field is associated to constant 
Peierls phases (i.e. exp(±i7r/2)), which further requires 
that the Rabi frequencies do not depend on the variables 
i, y and also obey the relation 2 ~ ^12 (^r equally 
^i.e, = i^l,e,y, using the definitions (28)). 

We stress that this concrete scheme leads to the SOI 
studied in the previous sections and that it should open 
a QSH gap in an atomic setup. In order to observe 
the QSH phases induced by such a synthetic SOI in a 
cold atom experiment, these non-trivial Peierls phases 
need to be engineered in a fermionic lattice. Again, 
one can consider the atom-chip proposal of Ref. [21]: 
different sets of wires, aligned along x and y, should 
trap the states Cj and gj alternatively in both directions, 
as illustrated in Fig. 8c-d. Additional "Raman wires" 
should then trigger RF transitions and couple the states, 
producing the induced-hopping and associated phases 
described above. Another possibility would be to apply 
the superlattice methods of Refs. [61, 76, 77] to the Lieb 
lattice. Once the gauge field is synthesized, this setup 
needs to be superimposed with a state-independent 
Licb lattice yielding the desired total Hamiltonian 

Finally, we note that the scheme presented in this Sec- 
tion could be simplified in order to reproduce the Abelian 
Haldane model on the Licb lattice [6, 17]. In this case, 
only two internal states g and e would be needed, instead 
of four. A realization of the Haldane model on the Lieb 
lattice would lead to integer quantum Hall states [17]. 



C. Detection 

The TB Hamiltonian and its long wavelength approx- 
imation are valid for a Lieb lattice populated by single- 
component fermionic atoms, e.g. ^^K or ^Li. In this 
case the atomic collisions are negligible at low temper- 
ature [36]. From the experimental point of view, time- 
of-flight imaging via light absorption [78] can be used 
in order to detect the presence of massless fermions. 
The harmonic trap potential V{r) — TOCj^r^/2 confin- 
ing the fermionic cold atom gas is ramped down slowly 
enough for the atoms to stay adiabatically in the low- 
est band while their quasi-momentum is approximatively 
conserved. Under these conditions, free fermions expand 
with ballistic motion and, from the measured absorption 
images, it is possible [79, 80] to reconstruct the initial 
reciprocal-space density profile of the trapped gas. Then, 
the local density approximation is typically well satisfied 
and the local chemical potential can be assumed to vary 
with the radial coordinate as fi{r) = /iq — V(r), where /io 
is the chemical potential at the center of the trap. For a 
system of cold atoms at temperature T, the atomic den- 
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FIG. 9: (Color online) Panel (a) Atomic density as a function 
of the chemical potential. Panel (b) Density of states as a 
function of the chemical potential. The full line corresponds 
to Aso ~ 0.25t and the dashed line to Aso — 0. 

sity in the bulk is uniquely determined by the chemical 
potential 

'^(^) = ^E //(l^'"'^)^"^- (29) 

Here Sq is the area of the first Brillouin zone of the Lieb 
lattice, and /(k, a,^) = [exp[(i?Q(k) — fi)/kBT] + 1]"^ 
is the Fermi distribution function, where Ea(k) is the 
energy spectrum of the Lieb lattice, cf. Eq. (1). Figure 9a 
shows the atomic density p for the bulk as a function 
of the chemical potential /i. The contribution from the 
highly degenerate topological band (la) manifests itself 
at jU = as a sharp jump in the atomic density. This 
feature is specific to the Lieb lattice. For small finite 
/i we note that p evolves proportional to /i^, which 
reflects the linear dispersion of massless fermions near 
the band center as well as particle-hole symmetry. On 
the contrary, for values of p, close to the maximum or 
minimum of the energy band, i.e. far away from the 
band center (where the long wavelength approximation 
can no longer be applied), p varies proportional to p. 
When we consider a finite value of the SOI, a finite gap 
appears in the energy spectrum. This coincides with 
the horizontal segments in Fig. 9a. Moreover, we can 
observe that the atomic density in proximity of the SOI 
gap is not anymore a quadratic function of the chemical 
potential but is behaving linearly. This is a signature of 
the mass term introduced by the finite SOI. In Fig. 9b 
we report the density of states (DOS) for two values 
of Aso- Note the gap opening for Aso 7^ and the 
existence of two robust van Hove singularities at p = ±2. 

Let us stress the great similarity between the den- 
sity and DOS obtained here for the Lieb lattice and 
those reported in Ref. [39] concerning the T3 lattice. 
However, these two lattices differ in their topological 
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order since the SOI opens a trivial insulating gap in 
the Tjj lattice [22]. This fundamental difference can be 
emphasized by computing the density profile for 
the realistic finite-size system: in this case, gapless edge 
states contribute to the Lieb lattice atomic density and 
slightly tilt the SO plateaus. Note that this important 
effect, which is the direct signature of the topological 
phase, could only be observed for sufficiently small 
lattices (in which case the number of edge-states is 
not totally negligible compared to the number of bulk 
states). Another difficulty is that the edge-states could 
easily be destroyed by the harmonic potential: they 
should therefore be stabilized by sharp boundaries [20] 
or by designed interfaces [21] . 

Finally, the methods for detecting topological proper- 
ties, such as quantized Hall conductivity [79] and chiral 
edge-states [20, 21, 40] have been discussed recently and 
could be easily generalized to our Lieb lattice setup. 

V. SUMMARY 

We have investigated the fermionic properties of a face- 
centered square (Lieb) lattice. This peculiar system is 
characterized by the presence of a single Dirac cone at 
the center of the first Brillouin zone and a flat energy 
band at half filling. In particular we focused on the mod- 
ification of this exquisite energy spectrum in the presence 
of an external magnetic field and a next-nearcst-neighbor 
spin-orbit interaction. In the former case, we have shown 
the opening of multiple gaps due to the occurrence of 
Landau levels, which leads to the formation of two Hof- 
stadter butterflies separated by the robust flat band at 
zero energy. We have characterized the topological na- 
ture of these gaps by investigating the IQHE. Inside the 
two Hofstadter butterflies, the Hall conductivity is quan- 
tized and each gap is characterized either by positive or 
negative values. Importantly, we find that the energy 
gaps separating the two butterflies have a trivial topo- 
logical nature as they are characterized by zero Hall con- 
ductivity. This fact is a direct consequence of the flat 
band's trivial order. In the case of the spin-orbit interac- 
tion we have shown the opening of two symmetric gaps 
around the flat band. These are characterized by a non- 



trivial 7j2 topological phase leading to the quantum spin 
Hall effect. We have demonstrated the existence of heli- 
cal edge states for an abstract cylindrical system, as well 
as for the more realistic finite lattice. 

We have further investigated the Lieb lattice present- 
ing a long wavelength approximation for the system 
Hamiltonian around the single Dirac cone. We have 
shown that the Hamiltonian can be expressed in a rel- 
ativistic form - similarly to the honeycomb lattice case - 
which is characterized by a set of pseudo-spin matrices of 
size 3x3. These matrices fulfill the commutation relation 
of an angular momentum and describe a spin-1 particle. 
Within the same approximation we have demonstrated 
that the spin-orbit interaction simply results in a mass 
term within the quasirelativistic Hamiltonian. Further- 
more, we have inspected the properties of the Landau 
levels. In addition to the dispersion relation ruled by the 
square root of the Landau level index - as in the case of 
the honeycomb lattice - we have shown that there are 
two competing Landau levels at zero energy, which are 
related to the lattice topology and to the symmetry class 
of the Hamiltonian operator. 

This system and its associated properties could be en- 
gineered using fermionic cold atoms placed in an op- 
tical lattice resembling the Lieb lattice topology. We 
have proposed a method for implementing Abelian and 
non-Abelian synthetic gauge fields in order to simulate 
the presence of an external magnetic field and a next- 
nearest-neighbor spin-orbit interaction term. We have 
further shown that these synthetic fields trigger the open- 
ing of energy gaps while preserving the robust fiat band 
at half-filling, properties which can be directly deduced 
from atomic density measurements. In particular, we 
emphasized that the Lieb lattice is very well suited to re- 
produce the intrinsic spin-orbit term introduced in Ref. 
[7], which in this case, can be simply decomposed into 
nearest-neighbour hoppings on a square sublattice. 
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